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Abstract 



The precise details of how myosin-V coordinates the biochemical reactions 
and mechanical motions of its two head elements to engineer effective proces- 
sive molecular motion along actin filaments remain unresolved. We compare 
a quantitative kinetic model of the myosin-V walk, consisting of five basic 
states augmented by two further states to allow for futile hydrolysis and de- 
tachments, with experimental results for run lengths, velocities, and dwell 
times and their dependence on bulk nucleotide concentrations and external 
loads in both directions. The model reveals how myosin-V can use the inter- 
nal strain in the molecule to synchronise the motion of the head elements. 
Estimates for the rate constants in the reaction cycle and the internal strain 
energy are obtained by a computational comparison scheme involving an 
extensive exploration of the large parameter space. This scheme exploits 
the fact that we have obtained analytic results for our reaction network, 
e.g. for the velocity but also the run length, diffusion constant and fraction 
of backward steps. The agreement with experiment is often reasonable but 
some open problems are highlighted, in particular the inability of such a 
general model to reproduce the reported dependence of run length on ADP. 
The novel way that our approach explores parameter space means that any 
confirmed discrepancies should give new insights into the reaction network 
model. 

Key words: Myosin, motor, processivity, mechanism, model, optimisa- 
tion. 
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Introduction 

A myosin protein is an ATPase which gains enzymatic activity by attaching 
to an actin filament (0, H, 0) ■ The myosin proteins use the chemical energy 
released in ATP-hydrolysis to create directed mechanical motion. More 
than 100 proteins have been identified as belonging to the myosin super- 
family and they are organised into some 18 subgroups (0, 0). A myosin is 
identified by a conserved 80 kDa motor domain and it is usually assumed 
that all the myosin motor proteins share the same biochemical reaction 
pathway when hydrolysing ATP The most studied of the myosins is 
the non-processive (muscle) myosin II whose main reaction path is found 
to follow the classical Lymn-Taylor scheme (0) describing the correlation of 
mechanical and chemical events. 

Myosin-V is a dimeric protein involved in the intracellular transport of a 
variety of cargos. The neck region of the two head elements is three times the 
length of the corresponding region of the myosin II heads (0). It is generally 
assumed that the long neck region acts as a lever arm and the size of the 
protein makes it possible for myosin-V to walk hand-over-hand (0, H, 0, Hoh 
following the helical repeat of actin. Myosin-V was the first molecular motor 
shown to be processive along actin filaments |~Llh, Previously established 
processive motors like kinesin and dynein use microtubules as the track for 
their processive directed motion, or in the case of RNA polymerase, DNA. 
Naturally, there has been much interest both in confirming the processive 
motion of myosin-V and in gaining insight into the details of the molecule's 
motion using recent developments in experimental techniques 0, 0, 0, Ed, 
H, E3, El, 0, EE El, EE E IH EE EE EI El, EI El El) Similarly, several 



theoretical models of the myosin- V walk with different levels of detail have 
been proposed HE EE EE EE El, EI El, EE El, El, El EE) • 



How myosin-V coordinates the biochemical reactions and mechanical 
motions of the two head elements of the protein to become an effective 
processive molecular motor is an open question. While there is a general 
agreement on the enzymatic reaction path for a single myosin head, the 
details of how myosin-V keeps the two heads' reaction cycles in phase is still 
unresolved. Likewise, there is still a lack of understanding of how external 
forces directly influence the kinetic mechanism of the walk and what role 
the internal strain in the molecule plays. As a result of the large number 
of measurements made on myosin-V in the last years, this is an appropriate 
time to make more detailed, quantitative models of the myosin-V walk, and 
to see how such models compare with what is found experimentally. 

The outline of the article is as follows: we first establish our model for 
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the processive walk of myosin-V and explain the model's underlying as- 
sumptions. In the results section we present analytic results for our, rather 
general, reaction network. These include the velocity, run length, diffusion 
constant and fraction of backward steps. These observable quantities are 
functions that depend on a number of parameters (mainly characteristic en- 
ergies for the various transitions between states). The unknown parameters 
are estimated by optimising the agreement of the model with a chosen rep- 
resentative set of experimental results. This approach, which is based on 
defining a cost-function, is outlined in the third section. We then present 
quantitative results for the optimised parameters and compare them with 
available experimental data and end with a discussion. Agreement with ex- 
periment is generally reasonable, though we are unable to explain the trend 
of the run length with varying ADP concentration. We discuss how our 
approach might be used to explore such possible discrepancies between ex- 
isting reaction network models and experiment, generating new insights into 
the myosin V stepping mechanism. 

Model 

The reaction network for our model is sketched in Fig. The main reaction 
cycle contains states 1 to 5. In state 4 the rear head starts its reaction cycle 
by releasing ADP, while the front head will not react before the whole cycle 
has been completed making it the new rear head. In state 5 the rear head is 
in the rigor state with no nucleotide attached, and the rear head will detach 
rapidly from the actin filament when it reacts with ATP from the bulk. 
The mechanical motion that moves the rear head into the front position is 
found experimentally to happen in two steps (J22L Eil ). from position 1 to 
position 2, and then from position 2 to position 3. In state 3 the front head 
is weakly attached to actin, so only one head is strongly attached in states 
1 to 3, making these the more vulnerable states for total detachment of the 
molecule from actin. When inorganic phosphate Pi is released from the lead 
head, going from state 3 to state 4, the myosin binds strongly to actin and 
the lead head makes the so-called "powerstroke" . Since the rear head is still 
attached, the lead head will not achieve the usual post-powerstroke angle 
with the actin filament, which will cause internal strain in the molecule. 
This is consistent with electron microscopy (EM) images (jjj, E3) which show 
strongly bound myosin-V to take a position similar to a telemark-skier's 
stance, indicating large internal strain in the molecule. 

The model given above for the main walk follows fairly closely the inter- 
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mediate states proposed by Rief et al. (Jl2|) which accounts for the most 
important experimental findings for myosin-V |l3h. Note that there is 
still disagreement in the literature as to which states are present in the 
main reaction cycle for the processive motion of myosin-V (see for instance 
Refs. Q El, E3, S Q). A natural consequence of the microscopic 
size of the motor protein is that both mechanical and chemical effects are 
important. We expect the external force, due for instance to the viscous 
drag of a cargo vesicle, to most strongly influence the mechanical steps in 
the cycle which involve translational motion. 

The main reaction path is reversible in the sense that one could, in 
principle, run the ATP hydrolysis reaction in reverse by pulling the motor 
backwards with an external force. (This would not be a very effective way 
to produce ATP since the motor is kinetically tuned to move efficiently in 
one direction only). It is important for a model which aims to investigate 
the mechanism of a motor protein to also include the effect of futile cy- 
cles and detachment rates. Futile cycles consume ATP without creating 
any net movement of the molecule while detachment limits the run length 
of the molecule. If, in an evolutionary sense, myosin-V were to tune the 
reaction rates of the motor domain it might well seek to maximise the for- 
ward velocity (requiring "weak binding"). However, it should probably also 
seek to minimise the impact of futile cycles and detachment rates (requir- 
ing "strong binding"). The kinetic parameters that arise might then be 
expected to correspond to a suitable compromise solution. We remark that 
the computational scheme for exploring parameter space that we will later 
describe could be used in the future to explore such evolutionary pressures 
and trends. 

There will in general be a large number of possible unfavourable path- 
ways, but we will assume that among these there is one dominant futile 
reaction cycle which takes place when the ADP in state 4 of the main cycle 
(Fig. ^) detaches from the front head before the rear head moving to state 
6. In the futile cycle the front head then reacts with a new ATP molecule 
so that the molecule returns back to state 2 without having created any 
net movement. Similarly, we will identify what we believe to be the dom- 
inant detachment rate. We make the conjecture that the molecule is most 
vulnerable to detachment when only one head is attached to actin, as in 
state 2 (due to the high internal strain in the molecule, state 1 is a short 
lived unstable state). We will assume that detachment from state 2 domi- 
nates over all other detachment rates. The mechanism of detachment will 
be the release of ADP followed by the binding of ATP to the rear head so 
that myosin-V detaches completely from the actin filament before strong 
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attachment is achieved by the front head (state 7, Fig. [IJ. 

The assumption we have made so far is that the main cycle is described 
well by the qualitative model of Rief et al. |l2h where in addition the mechan- 
ical motion of the rear head to the front is taken in two steps. Furthermore, 
we assume that there is one futile cycle that dominates over all others, and 
that this futile cycle involves states 2 and 4 in the main cycle. Lastly, we 
assume that the detachment rate at state 2 dominates all others. The com- 
plete model we have outlined (shown in Figs. n an dEJ) is a minimal realistic 
model for the walk of a molecular motor like myosin-V. 



Reaction rates 

The reaction rates between the different states are described by Arrhenius 
expressions: 

Wt = r- 1 e -(G!+AG i )/* B r ) (1) 

where r is the fundamental time scale of the reaction and /cb is the Boltz- 
mann constant. (We use the notation where Ui and Wi are the forward and 
backward rates, respectively, away from state i (jH3)). G\ is the energy bar- 
rier between state i and its neighbour state in the forward direction, while 
AGi is the energy difference between the two states. We set T = 298 K in 
this article. The total energy balance for the ATP hydrolysis is 

AG hyd = fc B rin( [A [ ^J + J2 ■ ( 2 ) 

where the nucleotide concentrations are made dimensionless by dividing by 
the concentrations at the reference states for AGi, [ATP] = [ADP]° = [Pj 1 = 
1 M. The standard free energy AG (0) = J2 AG; ~ 32.5 [&J/mol] @) which 
is close to 13 k^T, while AGhyd — 25 k^T at cellular conditions (J2J). 

For the main cycle, one ends up with a total of ten reaction rates. The 
first four rate constants in the main cycle are related to the mechanical 
movement (see Figs . ^ and 0) ; 

ui = r^ 1 , (3) 

yj 2 = r -l e -(£strain ~ /ex (cW ~ | /ex/M VfoT ^ (4) 

m = T -l e -(Gl + fcx(d D + ±f cx /k H ) + bE stlain )/k B T^ ^ 
™ 3 = r^e-^^W, (6) 

where Td is a hydrodynamic time-scale related to diffusion over one step- 
length and /ex is the component of the external force parallel to the actin 
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filament. The mechanical step is separated into a so-called working stroke 
of dw — 25 nm and a diffusional sub-step of dp ~ 11 nm (see Fig. I3J), as 
indicated by experimental findings (|22l Eih . giving a total step size d = 
dw + do ~ 36 nm Thus we neglect the (weak) external force 

dependence of the diffusional step size in u\ but capture the dominant effect 
of the force in, e.g. retarding the activated rate U2 and in distorting the 
metastable shape of the bound arm in state 2 (through the ~ /g X terms in 
the exponents of W2 and u^). The total internal strain in the molecule -Entrain 
defines an effective Hookeian spring constant, kn, related to the compliance 
in the motor and neck region of the myosin head: 

-^strain = — ' (J) 

where we have assumed that the molecule is fully strained in state 1. 

State 2 is in mechanical equilibrium. The position of the hinge/neck in 
state 2 is influenced by the magnitude of the external force, which gives rise 
to the energy term 1/2/ c 2 x //ch in w-i and u 2 (Eqs. @ and ©)• When 
moving from state 2, either back to state 1 or forward to state 3, the 
molecule increases its internal strain by -Entrain ° r ^strain = l/2knd D re- 
spectively (Fig. [3]). Using the latter expression together with Eq. (JJJ) gives 
b = (dr)/d\y) 2 — 0.2. Since -Entrain > fr-E strain there is a bias in the forward 
direction away from state 2, making the molecule a Brownian ratchet Ipjij ). 

Notice that the external force / ex is defined to be positive in the direction 
opposite to the movement of the myosin-V molecule, and that only the 
influence of the force parallel to the actin filament is taken into account. 
(The external force is of course a vector quantity, with a possible important 
influence on the walk arising from the normal component of the force l|35l 
Ei3).) The external force will accelerate or slow down the two reaction rates 
W2 and U2 away from state 2, depending on the sign of f ex . 

Experimentally it is found that there is one rate- limiting step in the walk 
of myosin- At large external force, f ex , one or more sub-steps 



that couple to the force, will become rate-limiting. A phenomenological way 
to model this is (fllh 

n = n + t 2 e /ra< W feBT . (8) 

Here r* is the (average) dwell time for one step, while t\ and r 2 are the dwell 
times of the force independent and dependent sub-steps respectively, and the 
external force couples to a distance d e Q. When fitting experimental data to 
this phenomenological equation, naively one would expect d c s ~ 36 nm 
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which is the average step length. Instead it was found that d e g is between 
10-15 nm. From our model this scale emerges quite naturally, since the 
external force couples to the diffusional search governed by rate u 2 which 
gives r* ~ 7~i + l/u2- From Eq. (J5J) we have 



u 2 



where A is some constant. Here dp is 11 nm while the correction term 
\ /cx/^H of the position of the hinge in state 2 is of the order of 8 nm when 
close to stall force. It is satisfying that this is consistent with the value of 
d c f[ cited above and indicates that the way we include the two substeps in 
our model is reasonable. (See Ref. for a discussion of the limitations of 
Eq. ©). 

The reaction rates for states with both heads attached to actin (see 
Fig. are given by 



u 3 = T -i e - G 3/ fcBT , (10) 

WA = [Pj] T" 1 e -( G 3+ AG 3-(l-^s«n)/feBT i (n) 

u 4 = r- 1 e-^ T , (12) 

w 5 = [ADP]r- 1 e-( G l + AG4 )/ fcBT , (13) 

u 5 = [ATP] t' 1 e- G V k * T , (14) 

wx = T -i e -^ + AG5 )/ fcBT , (15) 

where the nucleotide concentrations and [Pi] in Eqs. (|11|) - (|14[) are given as 
dimensionless quantities (see comments under Eq. (|2*jl). r is a microscopic 
time scale related to the characteristic oscillation frequency of the protein. 
The free parameters in the model are the activation energies G|, the energy 
differences AGj, and the strain energies. When an estimate is made for 
these free parameters, the resulting energy landscape (Fig. 0} gives directly 
all the predicted reaction rates, and the model's predictions for myosin-V's 
velocity and run lengths along the actin filament (see Appendix EJ). Note 
that t is not really an independent variable since it can be absorbed into 
the activation energies; 1/r exp(Gj) = 1/tq exp(G| +ln(ro/r)). 

As a measurement of the deviation away from equilibrium one can in- 
troduce the parameter 

5 

r = TT = e AG hyd /fc B T e -f cx d/k B T ^ / 16 x 

j=l 1 
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which gives the thermodynamic driving force for the molecular motor (J3 



371 ) . To fulfill detailed balance we have T = 1 at equilibrium 1)371). As 
expected from the energy balance the external force appears in T as f ex d', 
this is the total work done by the motor when completing one reaction 
cycle with step-length d. From Eq. 1)16)1 it is clear how the external force 
shifts the apparent equilibrium constant of the hydrolysis reaction . This 
coupling between the external force and the free energy of the hydrolysis 
of ATP, gives directly the thermodynamic upper bound on the stall-force 
/stall = AGhyd/rf (— 2.8 pN at cellular conditions). (The presence of futile 
cycles will influence the upper bound on the stall-force, but this correction 
is found to be insignificant for our (best) model). 
The reaction rates for the futile cycle are given by 

U 4 ,6 = Ui e - aE ^in/k B T ^ ( 17 ) 

Wq a = W5 e-^w/feT ^ ( 18 ) 

1*6,2 = U 5 , (19) 
U> 9 g = Q- E ztia.in/kBT g-/ 8 x (do+2 fex/ks)/k B T ^ ^0) 

Here Uij and Wij are the reaction rates from state i to state j. Using 
the fact that the two head elements are structurally identical, we make the 
assumption that the front head has to have the neck at a similar angle to 
that between the rear neck and the actin filament in step 4 of the main 
cycle to make it possible for ADP to be released or bind to the myosin 
motor domain (Eqs. 1)17)1 and 1)18)) ) • To achieve this, the front head has 
to overcome the strain energy in the molecule in state 4 which creates an 
energy barrier aE strilb i n for the ADP release reaction governed by rate 
(see Fig.|3J). In this way myosin-V synchronises the biochemical reactions of 
the two heads. 

To move from state 2 to state 6 (Fig. [^1, the molecule has to increase the 
internal strain by Strain without the energy from the hydrolysis reaction, 
(Eq. 1)20)1 ). so there is a very low probability for the futile cycle to run in 
reverse (effectively becoming a useful reaction path way). The reaction rates 
in the futile cycle follow the reaction rates for the main cycle except for the 
extra barrier caused by the strain in the molecule, so there is only one new 
parameter that appears in the the model. In principle aE str a,m could be 
estimated if the elastic moduli of the different parts of the molecule were 
known, by making a detailed structural model based on EM measurements 
jioT ) and crystal structures (0,1^) of myosin-V (see also Vilfan (pfflh). We 
will not attempt to do this here, but leave entrain as an undetermined en- 
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ergy barrier, an energy barrier used by myosin-V to synchronise the reaction 
cycles of the two head elements. 

The detachment from state 2 takes place if the ADP detaches from the 
rear head before the front head becomes weakly bound. The relevant reac- 
tion rates are 



The external force will increase the detachment rate of the single head (an 
effect we neglect when both heads are strongly attached, which should be 
a good approximation when considering the strain level in the molecule). 
Pulling experiments on SI give an apparent interaction distance of 5 = 
2.4 nm (|2ll ) for the external force (Eq. ()23j0 . 

An extended model where we would also consider the reattachment rate 
of motors is of course possible (|4o[ ) but not particularly relevant since the 
local bulk concentration of motors and of actin target sites is usually not 
well-controlled in an experiment. Our main focus in this model will be 
single-molecule experiments and their predictions. 

In our model we use the fact that the two head regions of the myosin- 
V protein are identical and have identical biochemical reaction paths. It 
might seem at first that this would be an obstacle to the processive motion 
of the protein, since the reactions of the two heads must be out of phase to 
ensure that at all times at least one of the two heads is strongly attached to 
the actin filament. However, evidence has been found in EM experiments 
that the intra-molecular strain affects the two bound heads asymmetrically 
ijfil Eiih. The model outlined above shows explicitly how this asymmetrical 
strain can be used by the two heads to coordinate their reaction cycles and 
to minimise the impact of the futile cycles and detachment rates. 

Parameters 

The fundamental time scales, r and can be estimated using the Stokes- 
Einstein relation, is related to the diffusion time of the whole head 
element (with a "diameter" of ~ 30 nm) over the step size dw ~ 25 nm, 
which gives t& ~ 10~ 5 s. Similarly, r is related to movement of a few 
nanometres, where the relevant length scale is the thickness of the head 
element, giving r ~ 10~ 8 s. Note that any difference between r or tj in the 
different steps will be absorbed into the activation energies G\ (see comment 



W7,2 = W5 , 



U2,7 = u 4 , 



u-j = M5 e 



\U\5/k B T 



(21) 
(22) 
(23) 
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under Eq. 115(1 so we are only interested in the order of magnitude of these 
time scales. Similarly, the activation energies will also be modified by the 
bulk pH and ionic strength, even though the influence of these solution 
properties are not included explicitly into the model. 

Estimates of some of the activation energies are available through chem- 
ical kinetic reaction rate measurements, but we will not try to guess the 
values of the activation energies but leave them all as free parameters. The 
model of the myosin-V walk outlined above is a fairly detailed model con- 
taining seven different states and 13 reaction rates, still the number of unde- 
termined parameters is relatively low. We have the four activation energies, 
Crf, three independent energy differences AG{ (where one energy difference 
will be dependent because of Eq. ©), and finally two terms connected to the 
strain level in the molecule, -Entrain and aE stTa ,i n . We are therefore left with 
nine undetermined parameters in our model of the myosin-V walk, which 
is a small number considering that the model is able to predict not only 
the walk's dependence on the bulk concentrations of ATP, ADP, and Pi, 
but also how the external force couples to the walk. The model is also de- 
tailed in its prediction of how the internal strain influences the walk and the 
level of strain in the molecule. Our model is therefore over-determined with 
respect to the available experimental data and provides many measurable 
predictions. 

Optimisation 

To search the parameter space computationally we define a cost function 
which quantifies the agreement between the model and the data. There is 
clearly some subjectivity involved in choosing the terms in the cost function: 
for example, one might be more inclined to include data that have been 
confirmed by several different research groups. We cannot investigate all 
possible permutations of our cost function since the optimisation approach 
is quite computationally demanding. However, the 17 experimental features 
that we identify, and encode in the cost function, are those that we expect 
a good model for the walk of myosin-V to reproduce. These include trends 
observed under variation of nucleotide bulk concentrations and the external 
force. Each term in the cost function was included to search for such a 
trend or to restrict the value of the energy jumps AGj between substeps. 
Likewise, we specify variances to indicate how large a deviation from each 
target value we deem acceptable. The cost function, A, given in AppendixlBl 
is constructed as a simple sum of 17 terms and each term gives an O(l) 
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contribution when within the accepted experimental error (as defined by our 
chosen variance). We believe that our cost function is a minimal encoding 
of the most important experimental trends, but it is possible to extend or 
modify our approach by including further or different data points if desired. 

We do not choose arbitrarily what values the nine undetermined pa- 
rameters should take but instead search the parameter space for favourable 
combinations using a technique based on simulated annealing (j4~0 ) of the 
cost function. To make the search effective, we first evaluate the cost func- 
tion at 50 million points in a parameter space of nine times 25 k^T. (This 
is the energy available in the hydrolysis reaction, Eq. ([%]). under cellular 
conditions. The internal energy .E'strain and energy differences AGi cannot 
be larger than the energy in the hydrolysis reaction, while there are no such 
limits in principle on the kinetic parameters and a-E^train)- To make sure 
these points are evenly spread in the parameter space we use the Sobol 
quasi-random sequence (|4ll . l42h . The fifty points with the lowest cost func- 
tion from the Sobol sampling were then passed to a simulated annealing 
routine, where the cost function is the "energy" term. The number of steps 
in the simulated annealing is chosen so that it is equal to the number of 
steps in a random walk over the average distance between the Sobol points. 

An attractive feature of estimating the free parameters in the model like 
this is that it has some similarities to the way the molecular motor has tuned 
the same parameters through evolution to achieve physiologically required 
velocities and run lengths under variable cellular conditions. 

Results 

After optimising the parameters for our model (Table it was found that 
there is one rate-limiting step, = 14.7 s _1 , corresponding to the largest 
value of G\ for i = 2, 5, which is the release of ADP from the rear head 
in state 4, Fig. ^ The ADP release rate is in quantitative agreement with 
kinetic measurements jlllll^. which estimates between 

10 s" 1 to 20 s _1 . From the model (Fig. ^) it can be seen why ADP release 
is the crucial reaction step for the myosin-V walk, since both the futile cycle 
and detachment depend on it. By slowing down the ADP release, myosin- 
V achieves a larger duty ratio for the myosin head, but more importantly, 
also reduces the flux around the futile cycle and similarly the detachment 
rate. Since the rate limiting step also determines the average velocity of 
the motor, there is a trade off when tuning the rate of 114. It is satisfying 
that the optimisation scheme of our model is able to reproduce the kinetic 
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tuning found experimentally. It is worth noting that no optimal solutions 
were found in other parts of the parameter space, i.e., no other possible 
combination of reaction rates could be found with either a different rate- 
limiting step or more than one rate- limiting step. 

The optimised "best fit" model parameters, as given in Tabled gi ye rise 
to U5 = 0.3 /iM s _1 for the attachment rate of ATP to the rear head in state 
5, a prediction that is somewhat slower than experimental estimates of u§ 
in the range from 0.6-1.5 (M s' 1 Q. The release of P; is found to be 
fast when the front head is attached weakly to actin (state 3, Fig. ^) with 
us = 3200 s , an order of mag nitude above the lower bound from kinetic 
measurements 113 > 250s -1 1)14 ) • 

It is a somewhat subtle point that the velocities and run lengths that 
our model predicts (see Appendix E]) are independent of one aspect of that 
model, specifically whether the mechanical motion between state 3 and state 
4 is assumed to happen before or during the release of P;. This is in spite of 
the fact that the forms of Eqs. (|1(J|> and (|11|) are sensitive to this difference: 
the (1 — b)E stTa [ n term would appear in the exponential in 113 and not in 

if the mechanical transition occurred before the phosphate release rather 
than during the phosphate release, as we assume here. The reason for this 
is that the optimization process fixes only the observable rates («3 and W4). 
This appears to suggest that there may be some ambiguity in the value of 
the energy G\ that appears in these rates, depending on whether the motion 
between state 3 and state 4 happens before or during the release of P;, i.e. 
whether the peak in the energy landscape occurs closer to state 4 or state 
3 respectively. We can nonetheless determine the form of the equations 
by considering the kinetic measurements for P; release. These are carried 
out using single myosin heads without internal strain. If, for the intact two- 
headed molecule, the (1— i>)-E s tram term was instead placed in the exponential 
in U3 (and the peak in the energy landscape was close to state 4) then the 
rate for a strain-free transition, such as would be expected to be the case for 
single headed molecules, would be larger by a factor of exp[(l — fo) -Entrain]) 
giving rise to unrealistic phosphate release rates that would be four orders of 
magnitude higher than 3200 s _1 . This provides quantitative evidence that 
the power stroke of the front lever arm takes place substantially after Pj is 
released, as is illustrated in Fig. 0] where the energy barrier between state 3 
and state 4 is close to state 3. 

We estimate the internal strain, -Entrain = 12.8 k-gT which gives a rigidity 
k = 120 k^T nm = 500 pN nm 2 for the lever arm l)43h. in agreement with 
other estimates in the literature 0). This gives a Young modulus of Y = 
0.6 [GPa] assuming the myosin neck has an effective radius of 1 nm, which 
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is comparable to what is found in similar proteins (0). The strain barrier, 
a-E'strain; preventing the futile cycle was found in the optimisation scheme 
to be 5.4 k-gT, which slows down the release of ADP from the front head 
(state 4, Fig. ~ 50-200 times compared to the rear head under changing 
nucleotide concentrations in the bulk. This is consistent with measurements 
reported by Rosenfeld and Sweeney (|24h. 

Among the more important properties of myosin-V are its velocity and 
processivity along the actin filament. We are interested in how these two 
properties are influenced by the bulk concentrations of nucleotides and the 
presence of an external force in our model. In single molecule experiments, 
all of these parameters can be controlled and monitored, giving rise to direct 
measurements of their influence on the myosin-V walk. 

The velocity of myosin-V has been measured at between 200-500 nm/s 
(jllih. which is consistent with the magnitude of the velocities predicted by 
the model. The velocity appears to follow a Michealis-Menten like form 
in which it becomes independent of ATP at high ATP concentration, while 
linearly dependent on ATP at low ATP concentrations (Fig. [5J). There 
is similarly found to be a strong dependence of the velocity on the ADP 
concentration. P; only has a measurable influence on the velocity (and 
the run length) at very larg e excess concentrations, similarly to what is 
found experimentally (|l3U 23h. When comparing the velocity dependence on 
concentration of ATP and ADP directly with experimental measurements 
|23l ) , we find reasonable quantitative agreement (Fig. [5J) . 

Our model reproduces the trend of velocity with increasing external force 
(see Fig. EJ) found experimentally in Ref. ([251) . Note that our model only 
considers the external force parallel to the motion of myosin-V along actin, 
while in an experiment optical loads are usually applied both along the axis 
of movement and perpendicular to it (jsilliih. This might explain some of the 
discrepancies observed at large external forces. When applying a negative 
external force (pulling in the forward direction), the velocity is not found 
to increase significantly as has also been found by |2tih. When the pulling 
force in the direction of the motion exceeds 2 pN, we find that the average 
velocity decreases, since the molecule tends, increasingly, to be pulled off 
the actin track thereby populating state 7 in the reaction cycle (Fig. ^) and 
reducing the forward motion. This reduction in the average velocity at large 
negative force was not observed by Clemen et al. (|26h. 

Related to the velocity of myosin-V is the dwell time, t^, defined by 
Eq. (|49|) in Appendix which is the average time it takes before the 



molecule takes a forward step ((ill . I27T ). Comparing the prediction of the 



model for the dependence dwell time on force with experimental results 
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(|lll . 1251 ) reveals fair agreement for high ADP and high ATP concentration 
and the correct trend, though poor quantitative agreement, for high ATP 
but low ADP concentration (Fig. [TJ. However, the model is not able to 
reproduce the finding (fill . Eih that the dwell time becomes independent of 
force at low ATP concentration. 

An important parameter for a processive motor is the duty ratio (@), r^, 
the average proportion of the time the head is strongly attached to actin. 
The duty ratio of myosin-V is found to be close to 90 % (|2J), while our model 
predicts even higher duty ratio at low forces (Fig.|HJ). As expected, the duty 
ratio is reduced at higher external force, since the mechanical motion from 
state 1 to state 3 is slowed down, and becomes rate-limiting around 1.6 pN 
at saturating ATP concentrations, which is consistent with experimental 
results of 1.5 pN HQ. 

Baker et al. (J23J) have measured run lengths of myosin-V as a function 
of ADP and ATP concentrations. When looking at the dependence of the 
run length on the ATP concentration, our model predicts a decrease in run 
length when increasing the ATP concentration (Fig. , which is consistent 
with experimental results (j23T). The run length is found from our model 
by finding the eigenvalue of the matrix of rate constants which gives the 
slowest relaxation time in the system (see Appendix EJ). Since there is one 
dominating eigenvalue, we get a single exponential decay in run length. 

We find a nonmonotonic dependency of the run length on the ADP 
concentration, as is also found experimentally. (The nonmonotonicity leads 
to the crossing over of the [ADP] = 1 mM and [ADP] = 1 /iM lines on 
Fig. at low ATP concentrations.) However, while Baker et al. (j^) finds 
a strong increase in the run length when decreasing the ADP concentration 
below [ADP] = 1 mM, our model predicts a decrease in the run length 
(Fig. EJ. It is not immediately obvious why our model does not reproduce 
this experimental finding. It is clear that the model neglects many possible 
futile cycles and detachment rates and maybe even other possible useful 
reaction cycles. For instance in Baker et al. (|23j) it is suggested that myosin- 
V needs two useful reaction cycles to be able to function under variable 
conditions. Even though our model only has one main reaction cycle it 
seems to be able to reproduce many of the experimental findings under very 
different conditions, casting some doubt on this earlier assertion. It does not 
seem at all obvious that it would be of physiological advantage for myosin-V 
to reduce its run length under increase of the ADP concentration. This may 
suggest that futile cycles not included in our model play a role. Finally, this 
puzzle indicates that the ADP dependence of the run length is worthy of 
further experimental investigation. 



A Kinetic Model Describing the Processivity of Myosin-V 



15 



The run length drops of exponentially when increasing the force as shown 
in Fig. llUI The run length is also found to be very influenced by pulling in the 
forward direction (negative force). This can be understood by considering 
Fig. |HJ The duty ratio increases with an negative external force since the 
mechanical movement from state 1 to state 3 is accelerated. This also makes 
it less likely that the molecule will end up in state 7 and detach (Fig. 
At high negative external forces though, the run length is decreased, since 
the larger force is pulling the molecule off the actin track, Eq. (|23j) . This 
non-monotonic behaviour in the run length as a function of force was not 
observed in a recent study (|26h which found the run length of myosin-V to 
be fairly insensitive to both positive and negative external force over a large 
range of values. Note that our model probably exaggerates the increase 
in attachment rate of the front head to actin when pulled strongly in the 
forward direction. For a strong force the position of state 2 (Fig. |3J) might 
change so much that the front head no longer is in the target zone slowing 
down the diffusion to the target site. If a correction for this had been 
included in our model we believe that the run length increase for a negative 
external force would be reduced. 

Using Eq. 1481 Fig. shows the fraction of backwards steps is found to 
be low until around 2 pN, consistent with experimental findings EH. As 
is shown in Fig. 111! the model predicts that the fraction of backwards steps 
are larger at low ATP concentrations, again in agreement with experimental 
results QEj. 

A quantitative measure of the stochastic deviations from uniform constant- 
speed motion fFig|12|) is given by the so-called randomness ratio ( 27 , 145I ) , 
p = 2D /Vd, where D is the dispersion given by Eq. 1501 The reciprocal of p 
gives a measure of the number of rate limiting steps, and for different ATP 
concentrations and external forces, it is found that the model gives only 
one rate limiting step. This is somewhat different than earlier theoretical 
predictions and experimental measurements of the randomness ratio 

can be used to differentiate between theoretical models. At high forces p 
diverges because of the vanishing velocity close to stall force. 

A biochemical reaction network is expected to be robust to small changes 
in the kinetic parameters (f-iil ). This robustness is needed to tackle both the 
natural changes that occur inside a cell during its lifetime and the fact 
that cellular biochemical reaction networks are highly interconnected, so a 
perturbation in one affects many others. A simple and effective robustness 
test is to see how the motion of the motors is affected by changes of ±5% in 
the different parameters. This also gives information on which parameters 
have the largest influence on for instance the velocity and run length of 
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the molecule in the model. Fig. ^3 shows that G\ and to a lesser degree 
-^strain have the largest influence on the run length, while AG4 is the only 
parameter which gives significant changes in the velocity when perturbed. 
The importance of G\ and AG4 in controlling velocity and run length is to 
be expected, since these are the parameters determining the rate constants 
«4 and Wi for ADP release from and recapture by the rear head between 
states 4 and 5, which we have found to be the rate-limiting step. 

Since the temperature enters explicitly in our equations, via the thermal 
energy scale k^T, it is natural to look at how the temperature influences 
the run length and the velocity. Some caution should be noted though, 
since it is found that some of the free energy terms themselves are known to 
be directly dependent on temperature (|2J). Fig. El is obtained under the 
assumption that none of the energy terms depend strongly on temperature. 



Discussion 

Many reasonable and well justified models of the myosin-V's walk exist in 
the literature, but the majority of the models are of a qualitative nature 
and often introduce unnecessary (and uncontrolled) approximations when 
employed to obtain quantitative predictions. Of quantitative models in 
the spirit of this work, one interesting early study for myosin-V is due to 
Kolomeisky and Fisher (j2?l ) who considered a two state model which took 
some account of step-size variations. It was found that there is a substep 
in the walk of myosin-V creating another possible reaction pathway, as also 



suggested in some experiments ()25l . l47l ). One shortcoming in our model, 



is that the model does not allow for the possibility that the geometry or 
mechanics of the walk can change in different regimes, for instance at very 
different bulk concentrations or external forces. It is still an open question 
whether motor proteins change behaviour in a dramatic way in different 
regimes, but some evidence suggests that stepping length is influenced by 
force (|22l. l25l. l2fil) . Some recent papers have tried to explicitly calcu- 

late the different strain energies for the different step lengths, which would 
be helpful to clarify the detailed mechanism of the motion of myosin-V. One 
problem in such calculations is that several angles and rigidities are not well 
known so that the number of free parameters becomes very large. 

There is necessarily some arbitrariness as to how many and which dis- 
tinct states are included in a model, since the concept is strictly only a useful 
approximation to the complex, fluctuating motion of the protein. Several 
sub-states (and thereby sub-steps in the main reaction cycle) which have 
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been identified experimentally or on thermodynamic grounds 
are not included in this model. This means that the sub-reaction rates be- 
tween some of the states have been combined into one effective rate. Where 
possible we have chosen not to excessively "coarse-grain" the state-space, 
e.g. by classifying many states together into fewer, more broadly defined 
states. Philosophically this seems wise if one is not certain a priori that 
such coarse graining of state space won't reduce the models precision or 
predictive power. We also feel that it is of value that the model takes into 
account all the steps where the protein reacts with smaller bulk molecules, 
since bulk concentration is something that can be controlled experimentally. 
Similarly, we wanted to separate out the different steps that are expected to 
couple to the external force. As a result, the proposed model can make di- 
rect predictions as to how the external force produced by an optical tweezer, 
for instance, should change the behaviour of the myosin-V walk. It is any- 
way important to remember that the complexity of a model such as ours is 
not really a function of the number of states but rather the number of pa- 
rameters, which remains small, the interpretation of the parameters, which 
remains physically clear, and the reaction network topology. 

The elastic strain in the molecule plays several roles in our model; most 
importantly, the strain in the molecule synchronises the chemical reactions 
of the two myosin heads. This is vital for making myosin-V an effective 
processive motor |22h . The synchronisation is caused by the slowing down 
of the ADP release from the lead head compared to the rear head in state 
4, since the internal strain makes it less likely for the lead head to have the 
optimal angle relative to the actin filament necessary for ADP release. A 
related benefit of this is that the slowing down of ADP release from the lead 
head minimises the impact of the futile cycle. Also, the strain biases the 
detached head to stay in the "target zone" (state 2, Figs. ^ and HJ) when 
doing the (11 nm) biased diffusional search for the next attachment site on 
the actin filament. This bias both increases the velocity of the motor, since 
less time is used in the diffusional search for the next target site on the actin, 
and increases the processivity by decreasing the time spent with only one 
head strongly attached to the actin. 

The model does not contain any direct dissipation of free energy AGdiss • 
Such dissipation by friction and heat loss caused by the non-equilibrium 
motion of the protein would reduce the total useful work that the molec- 
ular motor could do by reducing the effective free energy of the hydrolysis 
reaction AG c g = AGh y d — AGdiss- We have not tried to quantify the size 
of AGdiss but assumed that it is relatively small. It would of course be pos- 
sible to include AGdiss as an undetermined parameter even though we have 



A Kinetic Model Describing the Processivity of Myosin-V 



18 



chosen not to do this here. The model does take into account the energy 
dissipated by the futile hydrolysis of ATP by the futile cycle (Fig. ^) which 
does no work and creates no net movement, but since a.E s train is estimated 
to be relative large, the release rate of ADP from the front head in state 4 
(Fig. ^) is reduced two orders of magnitude compared to the release of the 
rear head. This indicates that myosin- V is a tightly coupled motor under 
all conditions and the futile cycle will not dissipate significant energy. 

The starti ng p oint for our model was the qualitative model proposed 
by Rief et al. (|12l ) and the observation that myosin-V has a so-called pow- 
erstroke movement and a diffusional search. There are several alternative 
mechanisms suggested in the literature, although they typically have a lot in 
common since there is general agreement on the single head reaction mech- 
anism. One influential alternative model has been proposed by De La Cruz 
et al. (fijl ) where the strong attachment of the front head to actin is trig- 
gered by the release of the rear head. It is also possible to envision there 
being several parallel reaction paths followed by the motor, where all these 
paths contribute significantly to the forward motion. It would be natural to 
assume that different paths dominate in different regimes, which would then 
provide a strategy for the motor to function well under varying conditions. 
One interesting example (|2^ ) of such a parallel reaction path mechanism, 
involved combining the Rief et al. and the De La Cruz et al. models into a 
version with two-paths. Taking into account the number of proposed mod- 
els in the literature and then the possibility of different combinations of the 
various models, the number of permutations of possible models is clearly 
very large. We argue that a more quantitative analysis of the competing 
models would often be useful, e.g. using a computational scheme similar to 
that presented here. 

We have presented a moderately detailed model of the walk of myosin-V, 
and compared the accuracy of the model with experimental measurements 
for velocities, run lengths, and dwell times at different nucleotide concentra- 
tions and external force. Predictions are also made for as yet unmeasured 
quantities, such as the internal strain in the molecule and the randomness 
ratio. The model also clearly shows how the internal strain can be used by 
myosin-V to coordinate its forward motion. A clear advantage of our model 
is that the physical significance of parameters in the model is transparent. 
This transparency makes the model a useful reference for comparison to 
future experiments and aids in the identification of elements of the model 
which are accurate and elements that need refinement. As our work and 
other recent quantitative studies have demonstrated, it is possible to give 
a clear analysis of very detailed models. Even when experimental data are 
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individually not conclusive, the large number of existing experimental mea- 
surements available should help to differentiate the models when they are 
analysed in such detail. In turn this will be useful in clarifying the underlying 
mechanism of the myosin- V walk, which is still not precisely understood. 

We remark that ours is the first attempt to explore carefully whether a 
full reaction cycle, with a nearly complete list of relevant parameters, is or 
is not able to reproduce the trends observed over a variety of experimental 
studies. There are several useful outcomes from our study. With the excep- 
tion of the dependence of run length on ADP concentration our model does 
appear to be broadly consistent with the experimental observations. We are 
still unable to explain even the qualitative nature of the run length depen- 
dence on ADP concentration. This is a puzzle, which might signify a flaw 
in the accepted models and which we flag for future experimental attention. 
In addition, we find that phosphate release must be coupled to a release of 
energy over much of the transition between states 3 and 4 if the model is to 
produce results that are consistent with known rates of Pi release. 
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A Analytical Expressions 

These expressions were used to compare with experimental data. 

The master equations for the probabilities, Pj, of finding the molecule 



in state j, can be deduced from Fig. [21 to be 

Pi = u 5 P B + w 2 P 2 -(u 1 +w 1 )P 1 , (24) 

P 2 = U\P\ + W 3 P 3 - (u 2 + W 2 + U 2 j + W 2fi )P 2 + Uq, 2 Pq + ^7,2^7 (25) 

P 3 = u 2 P 2 + w 4 P 4 -(u 3 + w 3 )P 3 , (26) 

Pa = u 3 P 3 + w 5 P 5 - (u4 + W4 + it 4)6 ) Pi + w 6)4: P 6 , (27) 

P 5 = u 4 P 4 + w 1 P 1 -(u 5 + w 5 )P 5 , (28) 

P& = 1*4,6 Pa + w 2 ,6 P2 - (ue,2 + ^6,4) P& , (29) 

P 7 = u 2J P 2 - (w 7 , 2 + u 7 )P 7 . (30) 
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In matrix notation we have 

P = MP, (31) 

where M is the 7x7 reaction rate matrix and the jth component of the vector 
P is Pj. 

The probability of finding the molecule in one of states 1 to 7 is not 
conserved under equation (|3 1|) since the motor drops off the actin track at 
a rate wj from state 7 (equation (jSOJ). In order to solve equation (|3*Tj) it is 
helpful to renormalise, eliminating drop-off, to get a probability-conserving 
equation. This is done by writing 

Pj = — e~ xt Pj (32) 

and choosing the constants (fj such that the renormalised probabilities Pj 
satisfy a conservative set of equations. It is possible to show that this can 
be done (see Ref. (0)) if 



rrT 



ip = -\ip, (33) 



and that then the vector P of renormalised probabilities satisfies the equa- 
tion 

P = MP, (34) 

where M is a renormalised reaction-rate matrix with uj = and renor- 
malised rate constants 

Uj = Uj and Wj = Wj — 1 , (35) 

<Pj <Pj 

for j = 1 ... 5, where the index is periodic with period 5, and 

ifj ifj 
= u i,j— an d = Wjj— , (36) 

¥?i (Pi 

for the rate constants not on the main reaction cycle. 

The slowest, dominant eigenvalue, — Ao, is negative for reasonable rate 
constants, while the faster eigenvalues can be complex giving rise to fast 
oscillations in the reaction network. 

At long times, P tends to the steady-state solution of equation (|3"I|) (jioh . 
found by solving M P = analytically. The relaxation times of oscillations 
are found to be fast compared to the stepping time and to 1/Ao, so it is clear 
from equation (|32|) that the solution for Pj at long times will be dominated 
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by the slowest eigenvalue, — Ao- We therefore choose cp to be the correspond- 
ing eigenvector of the transposed matrix M. T . Since there are seven states 
in our model, this eigenvalue problem is best solved numerically. 

In the following equations we assume all the rate constants are renor- 
malised and drop the tilde. We look for steady-state solutions Pj = 0. (For 
another approach to achieve analytical expressions for a fairly similar reac- 
tion network, see Kolomeisky (jHlh ). The net flux between two neighbouring 
states is given by J = PjUj — Pj + iWj + \, and we seek to express the steady- 
state probabilities Pj in terms of the futile flux Jf ut = PqUq : 2 — P2W2,& and 
the main reaction flux J^yd = V/d = P\u\ — PiWi- It can be shown that the 
steady-state solution is given by 



P, 



V 



for j = 1, ...,5, P 7 



r - i 

(^2,7/^7,2)^2, and 

r 2 



r j'7 + s i J fut 



r 2 -i 

where T is given by Eq. (|16|) and 



r{Jfat + s{ 



(37) 



(38) 



U^2 U2 u 3 n 4;6 
Wga w 2 ,e w 3 W4 



AG hyd /k B T 



3=1 3 

Here we have introduced the auxiliary functions 

fe 



u 



fc=li=l U 3 +l 



1 



U ; 



fc=l 



(39) 



(40) 
(41) 



for j = 1, 5, where 6j i is the Kronecker delta function (with 8ji = Sj+5,i) 
and where the indices on the rate constants are periodic with period 5. We 
have similar auxiliary functions relating to the futile cycle, namely 



1 



3 k w' ■ 

1 +1:113* 

fc=ii=i i+* 



W 3+i 



k=l 



i=l "i-H 



(42) 
(43) 
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for j = 1, 4, where w- = Wj and ttj = Uj except for w[ = w 2 



w' 



it6,2, and «4 



It! 



3 

1*4,6 



We now have a period of four, 



^2,6, 

i 



u 



'j+4' 



• +4 , and (5^ = <5j+ 4 ,i, since there are four states that participate in 
the futile cycle (Fig. QJ. 

Using the conservation of probability Yl Pj = 1 (since drop-off has been 
eliminated through renormalisation) and a second independent expression 
for P 4 , 



P 4 



To-1 



r 4 Jf ut + s 4 



r/ 



(44) 



one can derive the analytical expressions for the velocity, V, and the futile 
flux, J fut : 

V = Bdi 



Ji 



ut 



where 
1 

P 



rr 2 



j'=i r i 



D 



+ 



r(r 2 

r4 



1) r 2 (r-i; 



r 2 (r-i) r(r 2 -i; 



(45) 
(46) 



r 2 (r-i) r(r 2 -i 



S4 



+ 



E 



j=i s i 



+ 



r(r 2 - 1) 
n 



r 2 (r 



with f,- = rj and Sj 



r 2 (r-i) r(r 2 -i); Vr 2 (r-i) r(r 2 -i). 

(47) 

for j = 1,3,4,5, while r 2 = r 2 (l + u 2 j/wt^) and 



Si 



s 2 = s 2 {l + u 2 j/w 7 ^)- 

With the probability of finding the molecule in each state now deter- 
mined, we can use previously established analytical expressions. The frac- 



tion of backward steps is ((23, E3) 



7T_ 



WqW\ 



UqUi + WqW\ 



(48) 



where u = u 5 (P 5 / Yl]=2 P j) and w = w 2 (P 2 / 



The mean forward-step dwell time is given by 

u + u\ + w + wi 

td = ~ 

u ui + wqw± 

and the dispersion is given by (jj^ . \bc 



, 5. 



(49) 



D = l 
2 



WqWi 



-1-2 



1 



(u + U1+W0 + W])' 



W Q W X 



(u + U\ + Wo + Wl) 

(50) 
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Since the probability of the myosin remaining attached to the actin fil- 
ament decays exponentially with the dominant eigenvalue — Ao, the typical 
duration of a run is 1/Aq and hence the run length L is given by 



B Cost Function 



The cost function contains 17 terms. Except where other concentrations are 
explicitly mentioned, the nucleotide concentrations are given by [ATP] = 1 
mM and [ADP] =[P;] =0.1 /iM. Likewise, the external force is assumed to 
be zero, except if stated otherwise. 

The first term in the cost function is constructed from the model velocity, 
V, compared to the velocity V E estimated from experimental data in the 
literature (fill . IT3 . H-l . IIR EiL Eih . The mean squared uncertainty in the 
measured velocities, <Ty E , are in general found to be around 10%. 

a (l) _ [V-V E ] 2 _ [^- 540nm/s] 2 

a\ E (54nm/s) 2 " 1 ' 

When the ADP concentration increases the velocity V E is found to decrease 
UH 0, El We choose a second and third ADP concentration, 

[ADP]( 2 ) = 200 fjM and [ADP]® = 2.5 mM, and construct additional terms 
in the cost function at these conditions: 



(2) _ [T/([ADP]( 2 ))-T/ £ ([ADP]( 2 ))] 2 _ [FQADPp)) - 320nm/s] 2 

(53) 



a vE([ADP]W) (32nm/s) 2 



a(3) _ [F([ADP]( 3 )) - T/ £ ([ADP]( 3 ))] 2 _ [V([ADP]( 3 >) - 130nm/s] 2 
4^([ADP](3)) (13nm/s) 2 

(54) 

When the ATP concentration decreases the velocity V E is found to decrease 
,m). Choosing [ATP]® = 10 fjM we construct a fourth term: 

(4) _ [VQATP]®) - V E ([ATP]W)] 2 _ [V([ATP]W) - 75nm/s] 2 
CT ^([ATP](2)) (lOnm/s) 2 

(55) 



A Kinetic Model Describing the Processivity of Myosin-V 



24 



It is found that changing [P,! several millimolar does not significantly change 
the velocity of myosin- V QUI). ChoosingjPip) = 40 mM, we add a fifth 
term using measurements by Baker et al. (|23i): 

A ( 5 ) = [^([Pi] (2) )-^([Pi] (2) )] 2 = [y([Pi]( 2 ))- 440nm/s] 2 ^ 
<4s ([ p. ]( 2)) (44nm/s) 2 

Baker et al. (j23h have performed a large number of experiments of the run 
length L at different nucleotide concentrations finding that the run length 
is increased with a decrease in the concentration of either ADP or ATP. We 
introduce four terms to the cost function based on run lengths at different 
nucleotide concentrations (where [ATP]^ 3 ) = 100 uM): 

A (6) = [L~L E ] 2 [L-800nm] 2 

a\ E (150 nm) 2 ' 1 ' 



(7) _ [L([ADP]( 3 )) -L E ([ADP]( 3 ))] 2 _ [L([ADP]( 3 )) - 400nm] 2 

7 L E ([ADP](3)) 



^^FAnPiCs). (150 nm) 2 



(58) 



a(8) _ [L([ATP]( 3 )) -L B ([ATP]( 3 ))] 2 _ [LQATPp)) - 1150nm] 
°L«([ATP](3)) (150 nm) 2 

A(9) = [^([Pi] (2) )-^([Pi] (2) )] 2 = [L([P i ]( 2 ))- 500nm] 2 
°^([Pi](2)) (150 nm) 2 

(2) 

The velocit y is found to be independent of external force up to / cx = 0.75 



(59) 



(60) 



pN (|12l. Il3l. l25r> and we include an tenth term based on this: 

f(2) r A r,Pl(2h _ mE(A2) r inT) l[2)M2 UHfW 



A(10 ) = [ADP]< 2 >) - F^/g, [ADP]( 2 ))] Z = [y(/j 2) )-y([ADP]( 2 ))]- 

^-Cfg?) (50nm/s) 2 

(61) 

where 1/([ADP]( 2 )) is the velocity in Eq. (|53j) . Similarly, it was found that 
run length is fairly independent of external force ([26T). giving the eleventh 
term: 



[L(/e ( x 2) ,[ADP]( 2 ))-L s (/ e (2) ,[ADP]( 2 ))] 2 _ [L(/ C (2) )-M[ADP]( 2 ); 



cj 2 (150 nm) 2 



(62) 
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where L([ADP]( 2 )) = 400 nm (23). Also dwell time is found to be indepen- 
dent of external force up to /Jx' = 1 pN (jlll . Eih giving rise to the twelfth 
term: 

( fi. 3) . r a dpi ( 2 ) ^ - tM (\ a dpi ( 2 ) W 2 ( fi 3) ^ - n. 1 r *1 2 

(63) 



A(12) = [t d (/ix 3) ,[ADP]( 2 ))-tf([ADP](2))] 2 [t d (/ c ( x 3) )-0.15 S ] 2 
<d„<3K (0.1 ^) 2 



Rosenfeld and Sweeney (J2J) found that the release of ADP from the front 
head is at least 50 times slower than from the rear head, which gives rise to 
the thirteenth term: 

A (i3 ) = _L(M 2 (64) 

50 2 V v ' y ' 

where Jf u t is the futile flux (Eq. (|4f)|)). The last four terms are restrictions 
on the possible values of the energy jumps AGj, reflecting some inherent 
limits caused by strict limits on the energy available in each sub-step of the 
chemical reaction. Using measurements on SI (see tabled as rough 'target' 
values, but allowing for a deviation from these values of a Ad = 3 k^T the 
last four terms are: 



A( i4-i7) = - ^gn 



(65) 

a Ad 

The total cost function is simply defined to be the sum of all the different 
cost terms: 

17 

A = ^A (i) . (66) 

i=l 

All the terms in the cost function cannot be expected to constrain completely 
independent properties of the model. However, with 17 differing constraints 
on only nine free parameters it is encouraging, and perhaps not surprising, 
that the optimization reveals that the best solutions are grouped in the same 
region of parameter space. 
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G\ G\ G\ G\ E strain aE strain AG 2 AG 3 AG 4 AG 5 



0.3 10.4 15.7 5.8 12.8 5.4 0.14 9.9 -10 13.1 

Table 1: The estimated values (in units of k-gT ) of the free parameters in 
the model as given by the optimisation routine (see also Fig. 0J. 



AGf AGf AGf AGf 
2 5.7 -7.7 15.3 



Table 2: Energy jumps (in k-gT units) measured for the reaction of SI with 
actin (taken from Table 14.2, Ref. (0) using Eq.EJ). 
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Figure Legends 

Figure ^ 

Sketch of the complete reaction network of the model. The Y-shaped molecule 
is the myosin- V protein which walks on actin filaments. The black actin 
monomers indicate the attachment sites spaced at ~ 36 nm. The labels T, 
D, and P; stand for ATP, ADP, and inorganic phosphate respectively being 
bound to the head. 

Figure 

The main and futile cycles combined in one scheme showing all the reactions 
paths between the 7 states in the complete model. See also Fig. ^ The 
reaction rates are given by the corresponding equations in the text. Reaction 
steps which release and bind ADP, Pj, or ATP are indicated. 

Figure Q3 

The mechanical movement of myosin- V takes place in two separate steps. 
The first step, through a distance dw — 25 nm, is from the highly strained 
state 1 to state 2 where the internal strain balances the external force. When 
the molecule diffuses to state 3, through a further distance g?£> ~ 11 nm, the 
internal strain increases to 6-E s train- 

Figure |3J 

The one dimensional energy landscape that we find for the walk of myosin-V 
in which the states in the model are indicated by the filled circles. Energy is 
measured in units of k^T . The energy changes associated with the dashed 
transitions are -Entrain and bE stTll m, being the energy barriers involved in 
moving away from state 2. Also shown is the rate limiting activation en- 
ergy G\ between state 4 and state 5. The generalised reaction coordinate 
X can be thought of as measuring the progress around the main reaction 
cycle (Fig. ^) . As such it reflects a combination of physical motion and the 
progress of biochemical reactions, according to the substep. The shape of 
the curve is somewhat arbitrary, but the peaks and the troughs are at the 
correct energies determined by the optimal values (Gj, -Entrain, and AGj). 
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Figure 

Predictions of velocity of myosin-V as a function of ADP, ATP, and Pj 
concentration that arise from our optimised (best) model with parameter 
values as shown in Table 1 (and used in all subsequent figures). In each 
case the other two reference concentration are taken from [ATP] = 1 mM, 
[ADP] =0.1 [M, or [P;] = 0.1 aM. The experimental data for varying [ATP] 
(squares) and [ADP] (circles) are from Baker et al. l|23h. 

Figure 

The velocity as a function of force. The solid line (and circles) when [ATP] 
= 1 mM and [ADP] = 200 fjM., the dashed line (and squares) when [ATP] = 
1 mM and [ADP] = 1 fjM, while for the dotted line (and diamonds) [ATP] 
= 10 fjM. and [ADP] = 1 fiM. The model (the lines) show similar trends to 
what is found experimentally (the circles, squares, and diamonds) |25h . 

Figure [7J 

Dwell time for [ATP] = 2 mM (solid line/circles), [ATP] = 10 ^M (dotted 
line/squares), and [ATP] = 1 mM and [ADP] = 200 /iM (dashed line/triangles). 
The experimental data are from Mehta et al. (jTH ) (circles) and Uemura et al. 
Hi]) (squares and triangles). 

Figure [SJ 

The duty ratio, r<j, as a function of force for [ATP]= 1 mM, [ADP] = [P;]=0.1 
fiM. 

Figure Q5J 

Run length L for different concentrations of ATP when the ADP concen- 
tration is equal to 1 mM (dotted-dashed line), 100 (j,M (dashed line), and 
10 fiM (solid line). At low ADP concentration, the run length becomes in- 
dependent of ATP concentration. The prediction of the model is compared 
with experimental results at low ADP concentrations (circles) 

Figure I10L 

Run length L for different strengths of the external force when [ATP] = 1 
mM and [ADP] = 200 fiM. For negative external force the run length has 
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a non-monotonic behaviour, where it increases ten fold before decreasing 
again. 

Figure 1111 

The fraction of backward steps, Prev/pfon is insignificant in the model until 
an external force of ~2 pN is reached. The solid line is for [ATP] = 2 mM and 
[ADP] = 200 fiM, while the dashed line is for a reduced ATP concentration 
of 100 fiM. 

Figure [T^ 

The randomness ratio, p, as a function of force at different ATP concentra- 
tions. 

Figure Q]3 

The circle shows the unperturbed velocity and run length for [ATP]= 1 
mM, [ADP] = [P;] =0.1 /iM and zero external force. The squares show the 
influence on run length and velocity of changes of ±5% in each of the nine 
free parameters of the model (while keeping the other parameters fixed). 
The run length is very sensitive to changes in G\. Also perturbing in -Entrain 
gives quite large change in run length. Large variation in velocity was only 
observed when perturbing the parameter AG4. 

Figure 1141 

Temperature dependence of velocity and run length (when [ATP] = 1 mM 
and [ADP] = 200 fjM.). The model predicts an increase in velocity with 
temperature, but a decrease of the run length. The velocity is found to be 
more sensitive to changes in temperature than the run length. 
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Figure 2: 
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Figure 6: 



A Kinetic Model Describing the Processivity of Myosin-V 



41 




Figure 7: 
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Figure 8: 
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